Phonons in pristine and imperfect two-dimensional soft colloidal crystals 
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The vibrational modes of pristine and polycrystaUine monolayer colloidal crystals composed of 
thermosensitive microgel particles are measured using video microscopy and covariance matrix anal- 
ysis. At low frequencies, the Debye relation for two dimensional harmonic crystals is observed in 
both crystal types; at higher frequencies, van Hove singularities in the phonon density of states are 
significantly smeared out by experimental noise and measurement statistics. We introduce methods 
to correct for these errors, which can be applied to disordered systems as well as crystalline ones, 
and we show that error correction leads to considerably more pronounced van Hove singularities. 
Finally, quasi-localized low-frequency modes in polycrystaUine two-dimensional colloidal crystals 
are identified and demonstrated to correlate with structural defects such as dislocations, suggesting 
that quasi-localized low-frequency phonon modes may be used to identify local regions vulnerable 
to rearrangements in crystalline as well as amorphous solids. 



PACS numbers: 82.70.Dd, 61.72.-y, 63.20.kp 

Recently, video microscopy has been cleverly employed 
to extract information about the dynamical matrix of 
ordered [T] and disordered colloidal systems, using 
a covariance matrix analysis. This technical advance 
has opened a novel experimental link between thermal 
colloids and traditional atomic and molecular materi- 
als [51-fT^. Along these lines, one ubiquitous feature of 
atomic glasses, the so-called "boson peak" due to an 
excess number of vibrational modes at low frequency 
[l3l 114] , has been observed in disordered colloidal pack- 
ings [H m [5]. Furthermore, connections have been es- 
tablished between "soft spots" associated with quasi- 
localized low-frequency vibrational modes and localized 
particle rearrangements in disordered colloids [H |9l ITSl - 
flT] . reinforcing the possibility that such phenomena 
might exist in atomic and molecular glasses, as well. 

The present paper has two primary themes. First, we 
use nearly perfect crystals to gain a better understand- 
ing of errors in the video-microsopy-covariance approach 
for studying phonon spectra, and to develop and demon- 
strate error correction approaches which are applicable 
even to disordered systems. Second, we study imper- 
fect crystals in order to probe directly the effects of de- 
fects on phonon modes. We find that structural defects 
in the imperfect two-dimensional colloidal crystals are 
spatially correlated with quasi-localized low-frequency 
phonon modes. Thus, our experiments extend ideas 
about quasi-localized low-frequency modes and flow de- 
fects in colloidal glasses [51 HZ] to the realm of colloidal 
crystals and suggest that phonon properties can be used 
to identify crystal defects which participate in the mate- 
rial's response to mechanical stress. 

Specifically, we employ video microscopy and covari- 



ance matrix analysis to explore the phonons of various 
two-dimensional soft colloidal crystals. By studying two- 
dimensional crystals [TS] , we avoid complications [HI [20] 
encountered by previous experiments which analyzed 
two-dimensional image slices within three-dimensional 
colloidal crystals to derive phonon properties [H [TU] . 
Our work is also complementary to earlier experiments 
by Keim, et al. ^ which studied two-dimensional crys- 
tals composed of particles interacting via magnetic-dipole 
forces and found good quantitative agreement between 
the dispersion relation measured by the covariance ma- 
trix technique and theoretical expectation. By contrast, 
the present experiments measure the vibrational density 
of states, which is much more susceptible to error, as well 
as the dispersion relations for crystals of particles with 
short-range interactions. In addition, we identify a nar- 
row band of modes in imperfect crystals that is associated 
with crystal defects. 

The experiments employed poly(iV-isopropylacryl- 
amide) or PNIPAM microgel particles, whose diameters 
decrease with increasing temperature. Particle diame- 
ters are measured to be 1.4 /im at 24 °C by dynamic 
light scattering with a polydispersity of 5%. PNIPAM 
particles are loaded between two coverslips. Crystalline 
regions are formed as the suspension is sheared by cap- 
illary forces. The samples are then hermetically sealed 
using optical glue (Norland 65) and thermally cycled be- 
tween 28 °C and room temperature for at least 24 hours 
to anneal away small defects. Particle softness permits 
the spheres to pack densely with stable contacts and yet 
still exhibit measurable thermal motions. Before data 
acquisition, samples equilibrate for 4 hours on the mi- 
croscope stage. Bright-field microscopy images are ac- 
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FIG. 1: (color online) Phonon modes in 2D colloidal crystals, a. Accumulated mode number, N{ljj), as a function of frequency, 
for a nearly perfect crystal (blue squares) and an imperfect one (red square); tj^ is drawn for comparison (black line), be. 
Spatial distribution of a low-frequency mode for the (b) nearly perfect and (c) imperfect crystal; the direction and magnitude 
of polarization vectors are represented by the direction and size of the arrows. 



quired at 60 frames per second with a total number of 
frames of 40,000. An image shutter speed of 1/4000 sec- 
ond is used. Each image contains about 3000 particles 
in the field of view. The trajectory of each particle in 
the video was then extracted using standard particle- 
tracking techniques [21]. Crystal quaUty is character- 
ized by Fourier transformation of the microscopy images, 
and by spatial correlations of the bond orientational or- 
der parameter, VPg |22j . To obtain intrinsic vibrational 
modes of the colloidal crystal samples, we employ covari- 
ance matrix analysis. Specifically, a displacement vector 
u{t) that contains the displacement components of all 
particles from their equilibrium positions is extracted for 
each frame. A covariance matrix C is constructed with 
Cij = {uiUj)t, where i and j run through all particles and 
coordination directions. To quadratic order, the stiffness 
matrix K, which contains the effective spring constants 
between particles, is proportional to the inverse of C with 
K = k-QTC^^ . Thus, by measuring the relative displace- 
ment of particles, interactions between particle pairs can 
be extracted and dynamical matrix constructed, i.e.. 



m mk^T 

The dynamical matrix yields the eigenfrequencies and 
eigenvectors of the shadow system: the system of par- 
ticles with the exact same interactions and geometry as 
the colloidal particles, but without damping. This ap- 
proach of measuring phonons permits direct comparison 
to theoretical models, e.g., models that might be used to 
understand atomic and molecular crystalline solids. 

The covariance matrix approach assumes harmonic po- 
tential energy landscape near the equilibrium configu- 
ration of the system. This assumption can be readily 
tested by the potential energy distribution of the ob- 
tained modes [52]. In our experiments, we find that 
some of the lowest frequency modes, typically less than 5 
modes, do not satisfy the harmonic assumption, possibly 



due to under-sampling of the lowest energy basins [23j . 
or more likely due to small sample drifts. These non- 
harmonic modes are excluded from our analysis. 

At low frequencies, Debye scaling requires that the 
phonon density of states D{oj) scales with w'*"^, where d 
is the dimensionality of the system, or that N{uj), the 
number of modes below a certain frequency w, scales 
as uj'^ [23]. Debye scaling for two-dimensional crystals 
is clearly exhibited by both our perfect and imperfect 
monolayer colloidal crystals (see Fig. [I^) . For more than 
one decade, N{uj) follows a power law close to 2, as ex- 
pected. The lowest frequency modes typically exhibit 
wave-like features as shown in Figs. [T|d and c (more real 
space vector plots of low frequency modes can be found 
in supplementary material |22j). Thus we conclude that 
the Debye scaling observed in Fig. [l^ is due to wave-like 
"sound modes." 

At higher frequencies, deviations from Debye behavior 
are observed in the phonon density of states, D{uj), as 
shown in Fig. [2^ for the pristine crystal. For triangular 
two-dimensional crystals with harmonic interactions, the 
density of states has two van Hove singularities, one for 
longitudinal modes and one for transverse modes (solid 
lines in Fig. [2^) ; these singularities are expected to arise 
at the boundary of the first Brillouin zone. In contrast. 
Fig. [2^ shows that the experimentally measured D{uj) 
(black open circles) exhibits a smooth peak at an in- 
termediate frequency and a shallow shoulder at higher 
frequency. 

We identify the peak and shoulder as vestiges of van 
Hove singularities. Several factors may contribute to the 
rounding of van Hove singularities in a colloidal crystal. 
For example, particle polydispersity may break transla- 
tional symmetry for the largest wave vectors, where van 
Hove singularities appear. The statistics associated with 
the finite number of frames (i.e., finite number of tem- 
poral measurements), as well as uncertainties in locating 
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FIG. 2: (color online) a. Phonon density of states of two di- 
mensional colloidal crystals. DOS obtained from experimen- 
tal data (black circles), DOS after suppressing non-nearest 
neighbor interactions in the stiffness matrix K (blue trian- 
gles), and DOS after N/T extrapolation(red squares). Nu- 
merically generated DOS for the harmonic triangular lattice 
(with matched sound speeds) is plotted as a guide of eye (solid 
lines), b. Dispersion curves for longitudinal (open symbols) 
and transverse modes (filled symbols) along high-symmetry 
directions, including uncorrected experimental data (black 
circles), data after non-nearest neighbor noise suppression 
(blue triangles), data after N/T extrapolation (red squares). 
Theoretical expectations are plotted in dashed (longitudinal) 
and solid (transverse) lines with matched colors 



particle positions, can introduce noise into the covariance 
matrix and thus into its eigenvalues and eigenvectors. In 
the following, we discuss these effects, and we show how 
to recover some of the expected behavior by applying 
corrections to the experimental data. 

One cause of the reduced or diminished peaks in the 
DOS is noise that gives rise to apparent interactions be- 
tween non-nearest neighbor particles in the stiffness ma- 
trix. In experiment, the effective spring constants be- 
tween non-nearest neighbor particles are never zero; in- 
stead they take on small values fluctuating around zero. 
These apparent interactions between distant particles are 
unlikely to be physical, since PNIPAM particles interact 
through direct soft sphere contacts. To examine the ef- 
fects of the spurious random interactions between distant 
particles, we set the matrix elements between non-nearest 



neighbor particles in the stiffness matrix K to be zero 
with corresponding corrections to the diagonal terms as 
required by translational symmetry, and then we solve 
for new eigenvalues. The blue triangles in Fig. [2^ show 
that suppression of the effects of these spurious interac- 
tions in the dynamical matrix leads to recovery of peaks 
at the frequencies expected for van Hove singularities, 
while maintaining Debye scaling at low frequencies. 

The apparent interactions between non-nearest neigh- 
bor particles are likely caused by a combination of two 
sources of error: uncertainty in particle displacement de- 
termination and limited statistics |25j. We note that 
in the limit of perfect statistics (i.e., where the covari- 
ance matrix is calculated from an infinite number of time 
frames), measurement error modifies the effective inter- 
actions in a way that is mathematically similar to the 
effect of quantum fluctuations at high temperature, as 
described by Wigner-Kirkwood theory, and to integration 
errors in molecular dynamics simulations |22j . Thus, in 
the absence of statistical error, measurement error would 
give rise to further neighbor interactions but would not 
smooth out the van Hove singularities. 

We now consider the opposite case, where there is no 
measurement error but there is error associated with the 
quality of statistics used to calculate the covariance ma- 
trix. A key quantity is R — N/T, where N is the number 
of degrees of freedom and T is the number of independent 
time frames used. Random matrix theory suggests that 
the eigenvalue distribution should converge to its limiting 
R — Q values linearly with R We have shown that 
the eigenvalues also converge linearly for crystals [52] . In 
the absence of measurement error, one could then expect 
to obtain van Hove singularities in the density of states 
by extrapolating to the limiting R = Q frequencies. The 
comparison between the raw data (black circles) and ex- 
trapolated data (red squares) in Fig. shows that the 
corrections are larger at higher frequencies, as expected, 
and that extrapolation modifies the shoulder at the sec- 
ond van Hove singularity into a small peak. 

We also studied the effect of errors on the dispersion 
relation, which is obtained from the eigenmodes of the 
covariance matrix. For each eigenmode, spectral func- 
tions for longitudinal and transverse components are cal- 
culated for high-symmetry directions [57| [55], and the 
wave vector corresponding to the maximum of the spec- 
tral function, kmax{^) is extracted. The binned disper- 
sion curves (black symbols in Fig. ^p) largely follow the 
theoretical expectation, obtained by fitting to the low fre- 
quency part of the curve to obtain the longitudinal and 
transverse speeds of sound, as shown in Fig.[2|D. However, 
the measured dispersion relation has a stronger depen- 
dence on /c, especially for the longitudinal branch. When 
we suppressed further neighbor interactions, the result- 
ing dispersion relation (blue symbols in Fig. [2j3) agrees 
significantly better with the theoretical expectation. We 
also extrapolated to the limit of perfect statistics, where 
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FIG. 3: (color online) Low frequency modes in a colloidal crystal with defects, a. A snapshot of an imperfect PNIPAM crystal 
with a grain boundary in the middle of the field of view. b. Participation ratio for eigenmodes in crystal with defects, c. Color 
contour plots indicate polarization magnitudes for each particle, summed over the low frequency modes with participation ratio 
less than 0.2. Circles indicate "defect" particles identified by local structural parameters. 



R — N/T approaches zero. We find from simulations 
that, like the mode frequencies, the dispersion relation 
approaches its limiting R — Q value linearly in R. Ex- 
trapolation of the data (red symbols in Fig. also leads 
to excellent agreement with theory. Thus, the dispersion 
relation appears far less sensitive than the density of vi- 
brational states to both of the leading sources of error in 
the covariance matrix technique: measurement error in 
the positions of particles and limited statistics. 

From the data corrected by suppressing further neigh- 
bor interactions and by extrapolating R to zero, we 
obtain (2D) bulk moduli of 2.2 • 10"^ kg/s^ and 1.8 • 
10^^ kg/s^, and shear moduli of 6.0 • 10^^ kg/s^ and 
5.6 • 10~^ kg/s^, respectively. The measured shear mod- 
ulus is in line with (3D) bulk rheology measurements of 
PNIPAM suspensions [^ . 

Finally, we explored phonons in the imperfect two- 
dimensional colloidal crystal shown in Fig. [3^. Most 
of the low frequency modes are extended and wave-like 
(as in Fig. [TJ;), consistent with the observation of De- 
bye scaling in the accumulated number of modes, Nitj), 
shown in Fig. [T^. To take a closer look, we calculated 
the participation ratio, which measures the de- 

gree of spatial localization of a mode uj is defined as 
Pi^) = (I]i"ii|e^,j|2) /{NY,imf\e^^i\'^), where e^,^ is 
the polarization vector in mode w and is the mass, of 
particle i. Fig. [3]d shows that while most of the modes 
have participation ratios near 0.5, the value expected for 
a plane wave, a few of the low-frequency modes have a 
significantly smaller participation ratio. 

The low-participation-ratio low-frequency modes ob- 
served in Fig. |3jD are reminiscent of quasi-localized low- 
frequency modes observed in colloidal glasses [5] . In 
jammed packings, such modes have unusually low barri- 
ers to rearrangements |16j and have been used to iden- 
tify a population of soft spots that serve as flow defects 
when the packings are sheared [T7]. Experiments on the 
PNIPAM colloidal glass similarly found that rearrange- 
ments induced by shrinking all the particles occur at soft 



spots [S]. In crystalline systems, it is known that re- 
arrangements tend to occur at crystal defects, partic- 
ularly dislocations. Our observation of quasi-localized 
modes in imperfect crystals therefore raises the question 
of whether the modes are spatially concentrated near 
structural defects such as dislocations. 

In the colored contour map in Fig. [3];, we plot the 
spatial distribution of the quasi-localized low-frequency 
modes with a participation ratio less than 0.2, i.e., 
J2 i^)^^ where No,2 is the number of modes 

pr{u)<0.2 

with participation ratio below 0.2, i is the particle num- 
ber, and is the polarization vector of particle i in mode 
UJ. The white circles in Fig. |3] indicate structural defects 
in the crystal sample, identified by local structural pa- 
rameters. Here a particle is identified as a "defect" par- 
ticle when the number of its nearest neighbors is not 6, 
or the magnitude of the local bond orientational order 
parameter — j^^k""^ e&'^^jk jg j^gg than 0.95. 

The spatial correlation between quasi-localized low- 
frequency modes (soft modes in crystals) and structural 
defects in colloidal crystals is obvious in Fig. [3j;. In 
particular, such modes in crystals appear to single out 
structural defects susceptible to external perturbations 
such as dislocations or interstitials; we note that they 
are less effective at picking out vacancies, which are me- 
chanically more stable. The correlation between quasi- 
localized modes and structural defects is robust to vari- 
ation of the participation ratio cutoff between 0.1 and 
0.3 m. 

The observation that quasi localized modes are con- 
centrated in regions prone to rearrangements not only 
in disordered solids [51[T7], but also in crystalline ones, 
suggests that such modes may be a general identifier of 
flow defects. Indeed, these modes, together with those 
in the boson peak at somewhat higher frequencies [30j . 
could serve as a useful link for understanding systems 
spanning the entire gamut from the perfect crystal to 
the most highly disordered glass. 
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